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Abstract —The simple Kalina cycle system 34 (KCS34)- 
has been studied to perform energy cogeneration from the 
waste heat recovery (WHR) in preheater cement 
industries. The preheater available energy was 
considered from a 5000 tc/day cement production 
capacity. Thermodynamic and simplified exergoeconomic 
models were developed in the Engineering Equation 
Solver (EES) software. Several cycle thermodynamic 
parameters as ammonia-water mixture concentration and 
turbine operating pressure were wide-ranging in order to 
maximize the cycle thermal efficiency aiming to minimize 
the electricity generation cost. The temperature-entropy 
KCS34schematics were shown for different best results 
aiming to understand which set of parameters targets the 
maximum KCS34performance. The produced power, the 
thermal cycle efficiency, the exergetic efficiency and the 
exergoeconomic electricity specific cost were plotted for 
the different ranges of the independent parameters. The 
optimum results for a range specific investment price 
were presented. The main conclusions indicate that in the 
range of the studied parameters the turbine operating 
pressure caused a generated power variation greater 
than the ammonia-water mixture concentration in the 
KCS34performance. It was also possible to conclude that 
the KCS34is competitive with the existing electricity 
prices. In this case the KC proved to be applicable for 
WHR in the cement industry. 

Keywords—Kalina Cycle, Waste heat recovery, cement 
industry, CHP, cost. 

I. INTRODUCTION 

Currently, the industrial sector is a pioneer in energy 
consumption. In Brazil there is a total of 88 cement plants 
according to the Cement Industry Union [1], which fueled 
an apparent cement consumption of 353 kg/person/year in 
2013, showing an increase of 1.4% over the previous 
year. Brazil is one of the world's largest cement producers 
placing fifth in the rank [1]. This way, the cement 
industry is in a large scale, a great investment option for 
energy recovery. Consequently, new technologies applied 
to this sector is a point to be thought out, discussed and 
worked mainly aim for improvements in industrial 


processes, the possibility of using thermodynamic cycles, 
especially the Kalina cycle in the utilization of waste heat 
during the cement production process for the cogeneration 
of electric energy [2]. The KCS34has been studied to 
perform energy cogeneration from the waste heat 
recovery in cyclonic preheater of the cement industries 
with dry production process. The preheater available 
energy was considered from a daily capacity of 5000 ton 
of clink. 

II. LITERATURE REVIEW 

Since the 1980s, when it was patented, the Kalina cycle 
has been an efficiency promising. The expected 
conversion efficiency is as high as 45% [3]. An 
optimization procedure for heat recovery boilers in 
combined power noticed higher exergetic efficiency with 
fluid ammonia-water compared with the use of pure fluids 
[4]. A Kalina cycle for electricity generation from the 
exhaust gases of a gas turbine combined cycle have 
attested that the Kalina cycle was found to be 10-20 % 
more efficient than the Rankine [5]. Other paper indicates 
that Kalina cycle has 3% higher performance than 
Rankine cycle in Husavick plant [6].In order to compute 
and locate the irreversibilities in a Kalina cycle, the 
exergetic efficiency values reached 55% for an input of 
the turbine steam temperature of 525°C and an ammonia 
fraction of 75% in the working fluid [7]. The exergy 
analysis of a cogeneration plant formed by a Kalina cycle 
and a four-stage desalination plant was performed with 
the purpose of the cogeneration plant was to generate 
electricity with the simultaneous production of fresh 
water from geothermal energy [8]. A detailed description 
of the Kalina cycle terminology for waste heat recovery in 
the cement industry has been presented by several 
authors. The advantages of using ammonia-water mixture 
as working fluid and points out the major design 
challenges for the application of waste heat recovery in 
the cement industry, which focus on the design of heat 
exchangers that recovers energy from gases and the 
design of the components of distillation and condensation 
system [9, 10]. Different Rankine cycle configurations 
were comparedfor waste heat recovery in the cement 


www.iiaers.com 


Page | 190 





International Journal of Advanced Engineering Research and Science (IJAERS) [Vol-5, Issue-11, Nov-2018] 

https://dx.doi.Org/10.22161/iiaers.5.ll.27 ISSN: 2349-649S(P) \ 2456-1908(0) 


industry [11]. The authors also mention that if you want 
to improve or increase the generation rate per ton of clink 
should choice to use the Kalina cycle because with it has 
a higher waste heat recovery in the cement industry. A 
Kalina cycle was optimized for waste heat recovery 
application [12]. A computer code was developed in 
Matlab to simulate the cycle. In the ammonia 
concentration range, between 0.8 and 1.0, at the turbine 
inlet, it was observed that a higher temperature in the 
separator and turbine inlet, which leads to increased 
efficiency, and that the maximum efficiency was obtained 
with a concentration of 0.9. An overview about the state 
of the art power generation technology from thermal 
sources temperature range with non-aqueous fluids shows 
the advantage of the Kalina cycle due to the evaporation 
process in the boiler at a variable temperature [13]. 

III. METHODOLOGY 

The study of waste heat recovery in the cement industry 
has been developed in the steps mentioned below: 

- Thermal schematic definition for waste heat recovery. 
The focus of this work is the cement plants that have 
satellite clinker coolers. In these plants the only available 
source of heat is the exhaust gas from the cyclonic 
preheaters; 

- Thermal cycle modeling. The modeling included the 
mass, energy, entropy and exergy balances. The estimated 
calculation of the cost of electricity generated was 
included into the model; 

- Heat recovery optimization. Different process variables 
taking place in the thermal cycle were found to maximize 
the generation of electricity and/or minimize the cost of 
the generated power. 

3.1 Thermal schematic definition 

Considering the thermal source in focus, was selected a 
KCS34 for the waste heat recovery. The schematic of the 
cycle is shown in Fig 1. Table 1 shows the data of the 
preheater available energy that was considered from a 
daily capacity of 5000 ton of clink. These data are at state 
‘17’ in Fig. 1. 

The Kalina cycle can be considered simply to have a 
small number of components. In the state ‘T it is 
specified the steam at the turbine inlet conditions ‘1’, 
specifically ammonia concentration, the pressure and 
temperature of the working fluid. The fluid has a rich 
concentration of ammonia at the turbine inlet in state ‘1’. 
After expansion, in the state ‘2’, the fluid is mixed with a 
poor solution coming from the ammonia separator - EPS 
(8). For mixing was employed a valve (9) to equalize the 
pressure of the mixture at the turbine exhaust. The leaner 
mixture to ammonia in state ‘3’ can be condensed at a 
lower temperature and pressure, thus allowing greater 
power generation in the cycle. Before entering the 


condenser (11), the heat exchanges operating fluid in the 
low temperature regenerator - LTR (6). In the pump (7) 
the fluid is pressurized to the operating pressure at the 
turbine inlet, hr LTR and high temperature regenerator - 
HTR (5) the working fluid is preheated before it enters 
the waste heat recovery boiler. The boiler is made up of 
the economizer - ECON (4), the evaporator - EVAP (3) 
and the super heater - SA (5). After passing the working 
fluid through the EVAP and ECON, the fluid enters SEP 
where it reaches the desired concentration of ammonia in 
the saturated vapor state ‘11’. The operating temperature 
at the turbine inlet is achieved by passing the fluid 
through SA. 


Table 1: Available energy for WHR [14] 


Composition 

% Mol 

C0 2 

28.0 

n 2 

69.0 

O 2 

3.0 

Parameter 

Value 

Temperature (K) 

623.15 

Pressure (kPa) 

101.32 

Volumetric flow (m 3 /h) 

269,526 



Fig. 1: The KCS34 


3.2 Thermal cycle modelling 

The modeling of the thermal cycle was performed in the 
Engineering Equation Solver (EES) software. The 
modeling included the mass, exergy, entropy and exergy 
balances as well as the calculation of exergetic Fuels and 
exergetic Products and exergetic efficiencies for each 
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component and thermal cycle as a whole. In equations 1 
to 6 we present the model for the Kalina cycle as a whole. 
In these equations, the different states correspond to those 
shown in the thermal cycle of Fig. 1. The different 
quantities and terms in these equations are: Q in , heat 
added to the Kalina cycle, in kW, Q out heat rejected from 
the Kalina cycle, in kW, W, Net generated power in the 
Kalina cycle, in kW, m, mass flow, in kg/s; h, specific 
enthalpy, in kJ/kg; ex, specific exergy (considering 
physical plus chemical, as explained later), in kJ/kg F, 
exergetic fuel, in kW; P, exergetic product, in kW; q ex , 
exergetic efficiency and t], thermal efficiency. 

Qin = ■ (h 10 - hg) + m n ■ (h t - h lt ) 

Qout = ™ 4 ■ (h 4 - h 3 ) 

W = m 1 - (h 1 — h 2 ) — rh 5 ■ (h 6 — h 5 ) 
vi' 

r] = T~ 

v ent 

_ P _ _ W _ 

^ ex F m 17 -(ex 17 -ex 20 ) 

k* = — 

Vex 


For the purposes of the calculation, we used the 
calculation subroutines gas mixtures properties presented 
by [15]. With these subroutines the properties in the states 
‘IT to ‘20’ are calculated. Specific exergy for the 
working fluid (states ‘1’ to ‘14’) considers both the 
physical and the chemical portions according to equations 
7 to 9. 

ex = ex f + ex q 
ex f = h - h 0 - T 0 ■ (s - s 0 ) 
ex q = x ■ ex q _ NH3 + (l - x) ■ ex q _ H; , 0 
In the last equations we considered T 0 = 295 K e P 0 = 
101.32 kPa to compute the values of h 0 e s 0 , x is the 
ammonia fraction in the working fluid. We considered 
341,250 U/kmol, 3,120 kJ/kmol and 11,710 kJ/kmol for 
the specific chemical exergy of the ammonia, the liquid 
water and the vapor water, respectively [16]. For the 
states ‘15’ to ‘20’ only the physical exergy was 
considered. Additionally, for modeling the data presented 
in Table 2 were assumed. In this table the temperature 
differentials presented refer to the state of higher 
temperature in the thermal cycle with respect to the state 
of lowest temperature, so the values are positive. In the 
case of A T sat 9 it is referred to the subcooling degrees 
under the saturation temperature at the turbine operating 
pressure. This guarantees that the working fluid will 
always be subcooled at the economizer outlet. The value 
of A T s represents the temperature difference between 
stages ‘17’ and ‘1’. The other fixed values are the cooling 
water inlet pressure (P 15 ), the ammonia concentration at 
the turbine inlet (xfj, and the pump (r) B ) and turbine (rj T ) 


isentropic efficiencies. The values that were presented in 
a range were used for optimization. 


Table 2. Input data for modeling 


Parameter 

Unit 

Value 

*i 

- 

0.96 

Pis 

kPa 

250 

ATsat_9 

K 

5 

&T S i 5 

K 

5 

*T 3 j 

K 

5 

A7-4_ 16 

K 

10 

A T s 

K 

20 

ATj 5 m 

K 

8 (1) 

A7'l2_13 

K 

44.7^2) 

Rb 

- 

0.85;3) 

It 

- 

0' 85 (4) 

Parameter 

Unit 

Range 

Pi 

kPa 

5,700 - CT00 

Qw 

- 

0.66 - 0K2 

*10 

- 

0.869 - 0925 

ATi9 9 

K 

11 -25 


The estimation of the power generated cost is done with 
equation (10) and data from Table 3: 

Cg er = CR 1 (c inv ' ~ + c q&m ) 


where 
C = 

'“‘inv 


1^1985 p / C ) 0,6 

IC2013 ' inV - f ' VW 


FA= 


i-(i+i) n 

(l+i) n -l 


(7) 

( 8 ) 
(9) 


Table 3. Input data for cost calculation 


Parameter 

Unit 

Value 

Reference 

CR 

- 

2.34 

[17] 

HO 

h/year 

6000 


i 

% year 

5 


n 

year 

20 


1985 

- 

175 

[18] 

7^2013 

- 

100 

[18] 

c ref 

kW 

6000 

[19] 

r. 

^inv re f 

US$/MW 

973,000 

[19] 

Cq&m 

US$/MWh 

0.096 

[19] 


The financial values were actualized and referenced for 
December 2013. In equations 10 to 12 C ger is the cost of 
the generated power, R$/MWh, CR is the change rate, in 
R$/US$, C lnv is the investment cost in US$/MW, FA is 
the amortization factor, HO is the operation time per year, 
C 0&M is the operation and maintenance cost, in US$/kWh, 
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i is the interest rate, n is the life time, IC is the cost index 
for a given year, C inl? is the investment cost at the 

reference capacity in US$/MW, c is the calculated 
generation capacity of the Kalina cycle, in kW, e c re f is 
the reference generation capacity, in kW. 

3.3 Heat recovery otimization 

The optimal region was located using the optimization 
procedure by genetic algorithms available in the EES. 

This optimization method was chosen because it is 
mathematically robust and running error-free with large 
number of variables and constraints. In order not to incur 
large calculation time was set to 3 the number of 
individuals, in 4 the number of generations, and in 2 the 
mutation rate. Mathematically the optimization problem 
was formulated by the Eq. (13). The values of the 
restrictions are higher than the values in Tab. 2 in some 
cases. The latter constraints are practical or conceptual, in 
the case of ^20 > it is the minimum temperature of the 
gases in the exhaust to atmosphere and in the case of the 
generation entropy, <j geri , which requires that in each of 
the system components meet the second law of 
thermodynamics. 

' Maximize W 

subject to: 

5,700 kPa < P 1 < 8,000 kPa 
0.40 < q 10 < 0.90 
< 0.20 < x 10 < 0.85 

0.90 <x ± < 0.96 
s 2s >0.0 kj/kgK 
T 20 > 340 K 

,a ger . > 0.0 kW/Ki = 1,2,..., 11 

IV. RESULTS AND DISCUSION 

Note that in Tab. 3 for most parameters the values found 
are close to one of their extreme limits. In the cases of P lt 
x 10 , x ± , and A T 12 13 the values were close to the 
maximum. In the case of A T 19 9 and A T s the values were 
close to the minimum. In the case of q 10 the value was 
close to the middle of the range. Further simulations 
showed that the variation of A T 12 13 and q 10 do cause 
variation in the generated power and the efficiencies of 
the Kalina cycle, hi the case of q 10 a very little increase in 
the generate power was noted at 0.74. On the other hand, 
the value of the A T s must be in the minimum to obtain the 
maximum power. At the optimal condition the net 
generated power was 2,725 kW, the thermal efficiency 
was 0.218, the exergetic efficiency was 0.552, a very high 
value which is explained by the input source of exergy to 
the system has a low exergy potential to be a residual gas 
stream of low temperature. At this condition the cost of 
the generated power was 278,03 R$/MWh, corroborating 
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the breakthrough that results using KCS34 for generating 
electricity from waste heat in the cement industry. 


Table 4. Values for the optimal condition 


Parameter 

Unit 

Value 

Pi 

kPa 

7,941 

<7io 

- 

0.6667 

*10 

- 

0.849 

*i 

- 

0.960 

A7l 2 JL3 

K 

47.48 

A7\ 9 9 

K 

11.73 

A T s 

K 

23.5 


The pressure of the working fluid in the turbine inlet has a 
strong influence on the net generated power. As we can 
see in Fig. 2 (A), as increases the pressure at the inlet to 
the turbine rise the net generated power. The net power 
also increases with increasing concentration of ammonia 
in the working fluid in the evaporator. Note that for the 
maximum concentration of ammonia that was studied in 
this simulation below the pressure of 6211 kPa at the 
turbine inlet the system cannot operate because violates 
the thermodynamics’ laws. The increase in the net 
generated power by the pressure is explained by the 
increase of the steam enthalpy at the turbine inlet. With 
increasing ammonia concentration at the evaporator outlet 
the generated steam flow increases. This happens because 
with more ammonia in the mixture the working fluid 
becomes zeotropic causing a uniform temperature profile 
in the evaporator, which reduces the destmction of exergy 
and allows you to generate more steam Note that(Ki)h the 
combination of these two parameters the net power varied 
hundreds of kW values in the ranges studied. 

In Fig. 2 (B) it is noted that with increased pressure at the 
turbine inlet, the thermal efficiency and exergy efficiency 
also increase. This increase is due to the increased 
enthalpy of the working fluid at the turbine inlet with 
increasing pressure. Note further that the influence of 
pressure on the efficiency is enhanced causing variation 
in several percentage points. 

Reducing the temperature difference between the 
evaporator outlet gas (state ‘19’) and the working fluid at 
the evaporator inlet (state ‘9’), i.e. the AT 19 9 , increases 
the generation of net power in simple cycle Kalina to 
increase the generation of steam to the turbine. Note in 
Fig. 3 (A) that the largest power values were achieved 
with the highest concentration of ammonia in the working 
fluid, reinforcing what was already explained before. It is 
also observed that although varied in the same proportion, 
the effect of increased ammonia concentration exceeds 
the reduction of the title at the evaporator outlet (which 
should reduce the generated steam flow). The ammonia 
concentration is therefore one of the most influential 
variables in the generation of net power in cycle. The 
temperature difference between the evaporator gas output 
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and the working fluid in the evaporator inlet A7\ 9 9 has a 
strong effect in the calculation cycle of the investment 
cost. This variable is directly bonded to the surface of the 
evaporator heat exchange, and thus to the cost. Thus, 
operating in the minimum level of AT 19 9 , will not 
necessarily lead to lower cost of generation. 
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Fig.2 The influence of the turbine inlet pressure on the 
KCS34 performance 


The influence of the temperature difference between the 
evaporator outlet and the gas working fluid at the 
evaporator inlet (AT 19 9 ) on the thermal efficiency and the 
exergetic efficiency can be observed in Fig. 3 (B). The 
thermal efficiency remains unchanged with the variation 
of this parameter because it depends mainly on the 
enthalpy at the inlet and outlet of the boiler and turbine, 
which is not altered by the A T 19 9 . On the other hand, the 
exergetic efficiency becomes sensitive to this parameter 
(AT 19 9 ), which once changed, inpacts on exergy 
destruction in the evaporator, which is the largest in the 
cycle. 
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Fig.3. The influence of the o A T lg 9 on the KCS34 
performance 


The effect of varying the pressure at the inlet to the 
turbine on the generating cost can be seen in Fig. 4 (A). 
According to this figure the cost of electricity generated 
increases with pressure, which has a practical sense, 
because to adopt higher working pressures requires 
greater material thickness and higher cost. But the right 
explanation for this behavior is the cost modelthat was 
adopted. As it was said before, the value of the cost of the 
electricity generated is associated to the generation 
amount and the same has been observed which increases 
with the pressure at the turbine inlet. 

Figure 4 (B) shows the effect of temperature variation in 
the evaporator difference (AT 19 9 ) on the generation cost. 
The cost of electricity generated decreases with increasing 
temperature difference in evaporator following the same 
trend behavior power according to the adopted cost 
calculation model. The trend shows that the cost model is 
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assertive in order to reflect that a greater temperature 
difference in the evaporator will reduce the purchased 
cost with thisequipment and thus causinga lower total cost 
of the cycle and in the generated electricity. 
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Fig.4. The influence of the turbine inlet pressure and o 
Ar l9 9 on the KCS34 generation cost 


From the information shown in Fig. 4 the generated 
power range between 2,500 and 2,850 kW with 
generation costs values ranging from 262.00 and 286.00 
R $/MWh. The optimization performed with genetic 
algorithms finds a middle ground in the region with 
power values and costs of 2,725 kW and 278.03 
R$/MWh, respectively. These results allow us to state that 
depending on the context of this great region it is possible 
to opt for a solution that leads to a minimum generation 
cost, with less generation of power, or a solution leading 
to a generation of maximum electricity but with a greater 
generation cost. In any case, the cost of generated power 
is competitive in the Brazilian market. According to the 
Ministry of Mines and Energy [20] the average tariff 
value supply practiced by CompanhiaEnergetica de Minas 


Gerais- CEMIG for the industrial sector in May 2014 was 
296.38 R$/MWh. 


V. CONCLUSION 

The study about the waste heat recovery in the cement 
industry with KCS34 lets you express the following 
conclusions: 

- The amount of electricity generated can reach 2,725 kW, 
with a thermal efficiency of 0.218, an exergetic efficiency 
of 0.552 and a generation cost of 278.03 R$/MWh; 

- The pressure of the working fluid in the turbine inlet has 
a strong influence on the net power generated. With 
greater pressure at the turbine inlet, the generated power, 
the thermal efficiency and the exergetic efficiency 
increase; 

- The net power generated also increases with increasing 
concentration of ammonia in the working fluid in the 
evaporator. The ammonia concentration is therefore the 
most influential variables in the generation of net power 
in the cycle; 

- With the combination of the pressure values at the 
turbine inlet and the ammonia concentration in the 
evaporator that were studied, the net power ranged 
hundreds of kW; 

- The influence of pressure on the efficiency is enhanced, 
causing variation in various percentage points; 

- Reducing the temperature difference,AT 19 9 , increases 
net power generation in KCS34. On the other hand, the 
exergetic efficiency becomes sensitive to this parameter 
since it impacts on exergy destruction in the evaporator; 

- The cost of the generated electricity increases with 
pressure and decreases with increasing temperature 
difference in evaporator following the same trend of the 
generated power; 

- The cost model is assertive to reflect that a greater 
temperature difference in the evaporator will reduce the 
purchased cost with this equipment and thus causing a 
lower total cost of the cycle and in the generated 
electricity; 

- hi the great region calculated the generation of values 
may vary in the range between 2,500 and 2,850 kW with 
generation costs ranging between 262.00 and 286.00 
R$/MWh. This cost range is competitive in the Brazilian 
market due to the value of the average supply tariff 
applied by CEMIG, which for the industrial sector in May 
2014 was 296.38 R$/MWh. 
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